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Problem of Photochemical Equilibrium of Ozone 
in Planetary Atmospheres: Ozone Distribution 
in the Lower Atiposphere of Mars 


By 

Shardanand and Gerald Grams* 

National Aeronautics and Space Administration 
Wallops Island, Virginia 

SUMMARY 

Several investigators have applied techniques used for the study of the ozone 
distribution in the terrestrial atmosphere to the atmospheres of other planets. This 
paper shows that such techniques cannot be applied indiscriminately to calculate 
photochemical equilibrium distribution of ozone in planetary atmospheres. Limitations 
associated with some of the earlier treatments of photochemical equilibrium distri- 
butions of ozone in planetary atmospheres are discussed, and a technique having more 
universal application is presented. In addition, ozone concentration profiles for 
the Martian atmosphere based on the results of the Mariner 4 radio occultation 
experiment and the more recent results from Mariner 6 and Mariner 7 have been 
calculated using this approach. 


INTRODUCTION 

The basis for quantitative discussion of the photochemical production of 
atmospheric ozone is the assumption of photochemical equilibrium whereby the number 
of ozone molecules formed equals the number destroyed in unit time and unit volume. 
In most studies of equilibrium distributions of ozone, the analysis treats only the 
photochemical reactions which would occur in an atmosphere of atomic and molecular 
oxygen (0 and C^) and ozone (0^) in the presence of inert molecules (M) . For this 
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case, the basic photochemical reactions are: 


0 2 + hv + 0 + 0 


A < 2424 A 


( 1 ) 


0 + 0 2 + M -* 0 3 + M k 2 


( 2 ) 


0 + o 3 + o 2 + 0 2 


(3) 


0^ + hv -> 0 + 0 2 


o 

a 3 A £ 11,800 A 


(4) 


0 + 0 + M + 0 2 + M k 5 (5) 

where a ! s represent the probability of dissociation of molecules and k T s represent 
the reaction rate constant. The third body, M, conserves energy and momentum in 
the reaction. It is not necessary to include reactions involving excited molecules 
whose contributions to the equilibrium concentrations are negligibly small (ref. 1) . 
The method of Wulf and Deming (refs. 2 & 3; see also Craig, ref. 4) for calculating 
equilibrium amounts of ozone is commonly used for analyses of the ozone distribution 
in the terrestrial atmosphere. This method utilizes reactions (1) through (4) 
only, the inclusion of reaction (5) is regarded to have a negligible effect on the 
numerical calculations for the altitudes of interest in the terrestrial atmosphere. 


The model of ozone photochemistry developed by Wulf and Deming has been applied 
in several studies of the atmospheres of Mars and Venus. For example, some con- 
clusions about the stratospheric circulation and wind patterns of Mars were derived 
by Paetzold (ref. 5) from his calculated vertical distribution of ozone. However, 

Marmo et al . (ref. 6) took exception to the use of this method for the study of the 
atmosphere of another planet and included reaction (5) in their calculations of the 
ozone distribution in the Martian atmosphere. Their results demonstrated the 
importance of the direct association of oxygen atoms by three-body collisions for 
the Martian atmosphere. In addition to a sharp reduction in the total amount of 
ozone in the atmosphere, significant differences in the shape of the ozone concentration 
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profile were obtained. The ozone concentration profile calculated by Marino and his 
co-workers increased mono toni cal ly with decreasing altitude to a maximum at the 
surface of the planet, as opposed to the ozone-concentration maximum at about 40 km 
predicted by Paetzold's analysis. Because of the importance of ozone in planetary 
atmospheres, the system of photochemical reactions leading to the equilibrium 
concentrations of ozone has been re-examined to better comprehend the implications 
of this discrepancy. 


Photochemical Scheme for Calculations of the 
Vertical Distribution of Atmospheric Ozone 


First consideration is given to the scheme of Wulf and Deming involving 

reactions (1) through (4) which will henceforth be referred to as the 4-reaction 

system. Let n n , n 0 , n OJ and n be the number densities of 0, 0 o , CL, and M, 

1 Z J m z j 

respectively. Then, under conditions of photochemical equilibrium, the rate of 
change of the concentrations of the various compounds of oxygen will be given by: 

dn i 

= 2ct n 0 + a,n - k 0 n.n o n - k n n, = 0 (6) 

dt zz jo ZlZm o±o 


dn 2 

= -a 0 n 0 + a 0 n 0 - k 0 n-n 0 n + 2k 0 n n n Q = 0 

dt z z jo zizm oij 


(7) 


dn 3 

= -an + k n n - k n n = 0 . (8) 

dt J J ZlZm olo 
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These three equations are not independent since continuity of mass considerations 
require that 


dn^ dn^ dn^ 
d t~ + 2 dt~ + 3 dtT 


(9) 


By a linear combination of equations (6), (7), and (8), the following two 
independent equations are obtained: 


a 2 n 2 = k 3 n l n 3 


( 10 ) 


a 2 n 2 + a 3 n 3 = k 2 n l n 2 n m 


( 11 ) 


When n^ is eliminated from equations (10) and (11), we obtain the quadratic equation, 


( a 3 k 3 ) n 3 + ( a 2 k 3 n 2) n 3 


a 2 k 2 n m n 2 


= 0 


( 12 ) 


In most previous studies, n^ has been calculated from equation (12) as a function 
of altitude to give a vertical profile of ozone concentration, using an assumed 
vertical distribution for the equilibrium values of n^. In the terrestrial 
atmosphere, this procedure has been well established for calculating the photo- 
chemical equilibrium distribution of ozone in the stratosphere, where most of the 
atmospheric ozone is known to exist. However, limited knowledge of the composition 
of the atmospheres of other planets prevents the application of such a simplified 
scheme as the 4-reaction system to studies of the ozone distribution in those 
planetary atmospheres. The earth T s atmosphere, with oxygen as a major constituent, 
is decidedly different from the atmospheres of Mars and Venus, for example, which 
are now regarded to be composed mostly of carbon dioxide, containing only traces 
of oxygen. Thus, a technique which has been successful in studies of the observed 
temporal and spatial variations of the ozone distribution in our atmosphere should 
not be expected to apply with equal success to the atmospheres of other planets. 

One aspect of the 4-reaction system that is particularly disturbing for studies of 
the ozone distribution is when equations (10) and (11) are solved to eliminate either 
n^ or n^. Then, the resulting quadratic equation in n^, 


( k 2 k 3 n m) n l " ( a 2 k 3) n l - a 2 a 3 = ° 


(13) 
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is independent of local concentrations of the other oxygen compounds. This implies 
that n^, n^, and n^ are no longer in equilibrium, contradicting the assumption of 
equilibrium conditions that led to the derivation of equation (13) • This contra- 
diction would indicate that the 4-reaction system does not represent a physically 
stable condition. Then, the question arises: why has the 4-reaction system been 
successful in calculating photochemical equilibrium distributions in the Earth’s 
stratosphere? This equation predicts that the number of oxygen atoms is determined 
only by the reaction rate constants, the local number density of all atmospheric 
molecules, and the probabilities of dissociation of oxygen and ozone molecules. In 
special cases, as the Earth’s stratosphere, and ot^ can provide an indirect coupling 
with the columnar number densities of 0^ and 0^ above the given altitude which, 
in turn, would be approximately proportional to the local values of n^ and n^. 

The values of and are evaluated by integrating over wavelengths X expressions 
of the form 


a 2 ( z ) 


/ 


2424 A 
( A) $(z , A) dX 


(14) 


and 


a 


3 



11,800 A 

3 (X) $(z ,X) dX 


o 


(15) 


where a ^ and are, respectively, the absorption cross-sections of 0^ and 0^, and $ 
is the solar flux at the altitude z. For the case in which the only absorbing 
constituents of the atmosphere are 0^ and 0^, the solar flux at normal incidence 
can be evaluated from the expression 


- ♦.(X) 


exp 


a 


(z) dz 


z 


o 3 (X) fn. 


(z) dz 


z 


(16) 
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where $ is the solar flux at the "top" of the atmosphere. If the columnar number 
densities of 0 ^ and 0^ are large enough to affect the numerical value of the 
exponential factor in equation (16), as is true in the Earth’s stratosphere, an 
indirect coupling, as described above, is possible. However, for columnar densities 
of O 2 and 0^ that are smaller than in the Earth’s stratosphere, the coupling no 
longer exists and the contradiction is present. For example, the equilibrium number 
density of oxygen atoms calculated from equation (13) at 100 km in the terrestrial 
atmosphere would be approximately 2 orders of magnitude higher than n^, the actual 
number density of all atmospheric molecules. The contradiction will also be present 
for planetary atmospheres in which 0 ^ and 0^ are no longer the major absorbers of the 
solar radiation. In that case, equation (16) would be rewritten with extra terms in 
the exponential factor which may be much larger than the terms that are related to 
absorption by and 0^ so that the coupling might again be destroyed. Since the 
atmospheres of Mars and Venus, for example, have oxygen present only as a trace 
constituent, it would be expected that the 4-reaction system is unsuitable for 
calculating equilibrium ozone concentrations in the atmosphere of these planets. 

In order to assess the problem in more detail, another aspect of the 4-reaction 
system is considered. In this part of the analysis, an equation for conservation of 
mass is introduced, 


n i + 2n 2 + 3n 3 = 2 ^2 ’ (17) 

where n^ is a constant which may be regarded as an original number density of oxygen 
molecules which, under photochemical action, has formed oxygen atoms, oxygen molecules, 
and ozone molecules. By combining equations (10), (11), and (17), an equation for the 
local ozone concentration is obtained: 


<*2 ( 2 ^ 2 _ “P 

3 3c* 2 + 2k^n^ 


(18) 
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Solution of this equation for the ozone number density requires a knowledge of the 
number density of oxygen atoms. By eliminating n^ and n^ from equations (10), (11), 
and (17), the following is obtained: 


( k 2 k 3 n m ) n i " (“2 k 2 + 2k 2 k 3Vm)4 " ( a 2 a 3 " 2a 2 k 3" 2 ) n l 


2a 2 «3 n 2 


0 , 


(19) 


which, at first glance, would appear to be coupled to the number of oxygen atoms 
available for establishing a condition of photochemical equilibrium. However, this 
equation can be factored in the form 


P k 2 k 3 n m' )n l “ ^ a 2 k 3^ n l 


a 2 a 3| 


3 ] [ 2n 2 n l] 


( 20 ) 


This result again contradicts the assumption of photochemical equilibrium. A 
solution based on the first factor is identical with the solution for equation (13) 
which has already been determined to be contrary to the physical meaning of equilibrium 
concentrations of oxygen compounds. The solution based on the second factor, 
n^ = 2 n 2 , represents complete dissociation of the available oxygen compounds which, 
again, contradicts the assumption of photochemical equilibrium. Thus, a conservation 
of mass relation cannot be incorporated meaningfully into the solution of the 4- 
reaction system. In view of the foregoing discussion, it has been concluded that the 
4-reaction system should not be applied without reservation for calculations of the 
ozone distribution of other planets and that investigations based on such calculations 
(see, for example, Paetzold (ref. 3) who studied the Martian atmosphere and Fabian 
and Libby (ref. 7) who studied the Cytherean atmosphere) should be reanalyzed within 
these limitations. 
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The system consisting of reactions (1) through (5) and, henceforth, referred to 
as the 5-reaction system is now considered. As before, the solution of the equations 
for the rate of change of n^, n^, and for equilibrium conditions leads to two 
independent equations. 


a 2 n 2 k 3 n l n 3 + k 5 n i n m 


and 


( 21 ) 


a 2 n 2 + a 3 n 3 


= VlVm 


+ k n n 
5 1m 


( 22 ) 


Equations (21) and (22), along with the conservation of mass relation represented by 
equation (17), are then available for solving for equilibrium number densities of the 
oxygen compounds. These equations may be combined to give the following expression 
for the number density of ozone: 


k 2 k 5 n l n m 


o 

a 9 a a + a 9 k 0 n - k 0 k 0 n^n 
2 3 ^ 3 1 23lm 


(23) 


If n^ and n^ are eliminated from the equations describing the 5-reaction system, a 
cubic equation in n^ results: 


^ k 2 k 3 n m “ 3k 2 k 5 n m " 2k 3 k 5 n m^ n l " ( a 2 k 3 + 2k 2 k 3 n m n 2 + 2a 3 k 5 n m^ n 
- ( a 2 a 3 " 2a 2 k 3^-) n l 


(24) 


+ 2a 2 a 3 n 2 = 0 
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The solution of equation (24) is coupled to the available oxygen as well as the other 
photochemical parameters and, therefore, provides the coupling mechanism for the 
desired condition of equilibrium. Thus, the 5-reaction system is to be preferred 
since, unlike the 4-reaction system, it does not contain an implicit contradiction of 
the assumed state of equilibrium. Furthermore, inclusion of reaction (5) has the 
advantage of increasing the generality of the technique for use in exploratory in- 
vestigations of photochemical equilibrium distributions of ozone in the atmosphere of 
other planets. Of course, the 5-reaction system is still not as general as a 
comprehensive analysis of all possible reactions such as the study of the photochemistry 
of a CO ^ atmosphere conducted by McElroy and Hunten (ref. 8) . However it is pointed 
out that, for the lower atmosphere of Mars from the surface to about 60 km, McElroy 

and Hunten assumed, as has been assumed here, that the relative concentration of 0 9 

-3 1 

is constant with a mixing ratio of 10 . While the 5-reaction system does not provide 

the amount of insight as a comprehensive study of the mechanisms which maintain this 
mixing ratio, it does not have intrinsic value for calculations of the ozone distri- 
bution in the lower atmosphere of Mars; and, therefore, such calculations have been 
performed based on the results of the Mariner 4 occultation experiment as well as 
the more recent results from Mariners 6 and 7. 


Calculations of the Ozone Distribution 
in the Atmosphere of Mars 


Calculations of the vertical distribution of ozone in the atmosphere of Mars 
based on the 4-reaction system and on the 5-reaction system have been reported by 
Paetzold (ref. 5) and Marino et al . (ref. 6), respectively. On the basis of the 
previous discussion, the results of Marino and his collaborators are more reliable. 
However, their atmospheric model has a surface pressure of 25 mb and partial pressures 
of 470 m STP for nitrogen, 55 m STP for carbon dioxide, and 70 cm STP for oxygen in 
accordance with the results of Kaplan et al . (ref. 9) and a temperature profile with 
a surface temperature of 250°K, a lapse rate of 2.9°K per km, and isothermal 
conditions above 30 km in accordance with the results of Schilling (ref. 10) . This 
atmospheric model is now outdated; furthermore, solar radiation was limited to the 
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spectral region 1750 to 3000 A in their calculations. For these reasons, it is now 
possible to provide better data of the ozone distribution of the atmosphere of Mars. 
Using an atmospheric model based on the results of Mariner 4 radio occultation 
experiment (ref. 11) and a computer program incorporating an iteration procedure to 
obtain self-consistent, steady-state distributions of ozone based on the 5-reaction 
system, calculations of the ozone content of the Martian atmosphere have been 
performed. In these calculations, no distinction has been made between 0(D) and 

O 1 

0( P) atoms as 0( D) are not present below about 50 km (ref. 8), the altitude interest 

of our major concern. The numerical data summarized in Table I were used for the 

calculations. The composition of the atmosphere was taken to be 99.9 percent carbon 

dioxide and 0.1 percent oxygen. Vertical profiles of 0, 0^ 9 and 0^ are shown in 

figure 1. Solid lines indicate calculated profiles based on the 5-reaction system 

and the dashed lines indicate profiles based on the 4-reaction system. The profiles 

for 0 graphically portray the problem regarding the 4-reaction system. For the 

5-reaction system, 0 profiles increase with altitude to the point at about 65 km at 

which almost complete dissociation of 0^ has occurred (n^s^n^). However, the 0 

profile for the 4-reaction system increases at a much faster rate and exceeds the 

value of n , which is taken to be the number density of all atmospheric molecules 
m 

at altitudes somewhat above ~60 km. For the 5-reaction system, the 0^ number density 
decreases monotonically with altitude and has an entirely different profile than the 
4-reaction system. Thus, the maximum concentration of Martian ozone is expected at 
the surface of the planet. Although the significant differences in atmospheric models 
results in a surface number density of ozone that is over an order of magnitude smaller 
than that calculated by Marmo et al . (ref. 6), it is concurred that an ozone maximum 
exists at the surface of the planet. This is in disagreement with the results of 
Paetzold (ref. 5) which indicated an ozone maximum at an altitude of about 40 km, 
based on the 4-reaction system and an atmospheric model similar to that used by Marmo 
and his co-workers. 

Additional calculations were performed to determine the effect of varying the 
parameters used in the model. None of the parameter variations yielded maximum 
concentrations of ozone at some altitude above the surface of the planet. The effect 
of lowering the content of carbon dioxide to 80 percent with 19.9 percent argon or 
nitrogen and 0.1 percent oxygen was to increase the ozone concentrations at all 
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altitudes by about 20 percent without appreciably altering the shape of the profile. 
When the oxygen content was increased to 1 percent with 99 percent carbon dioxide, 
the ozone concentrations were increased by a factor of about 2; a similar increase 
was found for an atmosphere of 80 percent carbon dioxide, 19 percent argon or 
nitrogen, and 1 percent oxygen. In addition, the effect of limiting the spectral 

o 

flux to the 1750-3000 A wavelength interval to match the calculations of Marmo 
et al . (ref. 6) was investigated. The resulting numerical values for the ozone 
number densities were about twice those obtained by using the solar flux from 

o 

1350-7100 A. This discrepancy can be attributed mainly to the dissociation of ozone 
by radiation in the visible part of the solar spectrum. 

Finally, atmospheric parameters obtained from Mariner 6 and 7 radio occult ation 
data have been utilized for calculations of the ozone distribution in the Martian 
atmosphere. The results of these calculations are displayed in figure 2 as vertical 
profiles of 0, 0^, and 0^ from the surface of Mars to an altitude of 35 km. The 
atmospheric parameters used in the calculations are those presented by Rasool et al . 
(ref. 12) for the Mariner 6 and Mariner 7 entry data which correspond to daytime 
conditions at low and high latitudes, respectively. Results obtained from Mariner 4 
data are also displayed for comparison; it is noted that the Mariner 6 and 7 results 
are not substantially different, with the exception of the calculated ozone concen- 
trations near the surface. Integrated total ozone amounts are calculated to be 
-3 -4 -4 

1.4 x 10 , 7.7 x 10 , and 8.8 x 10 cm STP for atmospheres corresponding, 

respectively, to Mariner 4, 6, and 7 results. It is pointed out, partially to give 
substance to the neglect of photochemical reactions involving CO^ and CO for cal- 
culations of equilibrium concentrations of ozone, that the calculated total ozone 
amounts are in excellent agreement with observations of the total ozone in the 
Martian atmosphere. The Mariner 6 and 7 ultraviolet spectrometer experiment detected 
ozone in the Martian atmosphere (ref. 13). Data from these measurements indicated 
a total ozone amount that was also about 3 orders of magnitude less than found on 
earth. These calculations, therefore, corroborate the measurements of ozone in the 
Martian atmosphere. 
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Hill I llll III! Ill III II 


I 


I 


III llll I llll I 


TABLE I. - PARAMETERS USED AND THEIR SOURCES 


PARAMETERS 


Mariner 4 


Surface temperature 


(a) All C0 2 — 180°K 

(b) 80% C0 2 and 20% N 2 and/or Ar — 175°K 
Surface number density 


(a) 

All 

C° 2 

— 1.9 x 10 17 cm - 

(b) 

80 % 

CM 

o 

O 

and 20% and/or 


Ar - 

- 2 

i -m 17 " 3 

. 1 x 10 cm 


Surface scale height — 9.0 km 
Temperature and density profiles 
Mariner 6 

Surface temperature — 250°K 
Surface number density — 1.7 x 10"^ cm 
Surface scale height — 12.5 km 
Temperature and density profiles 


SOURCES 


Fjeldbo et al . 
(ref. 11) 


Rasool et al . 
(ref. 12) 
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TABLE I. - PARAMETERS USED AND THEIR SOURCES - Continued 


PARAMETERS 

MARINER 7 

Surface temperature — 224°K 
Surface number density — 1.6 x 10^ cm 
Surface scale height — 11.2 km 
Temperature and density profiles 
Oxygen mixing ratio — 10 


Solar flux on the top of Earth’s atmosphere 
in 1350-7100 A region 

7100 - 3000 A 

3000 - 2500 A 

2500 - 1500 A 


1500 - 1350 A 

Average dilution factor to account for 
diminution of solar flux in vicinity 
of Mars — y = 0.4307 


SOURCES 


Rasool et al. 
(ref. 12) 


Belton & Hunten 
(ref. 14) 

McElroy & Hunten 
(ref. 8) 


Craig (ref. 4) 

Johnson (ref. 15) 

Detwiler et. al . 
(ref. 16) 

Watanabe (ref. 17) 
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TABLE I. - PARAMETERS USED AND THEIR SOURCES - Continued 


PARAMETERS 


SOURCES 


Absorption cross sections: 
0 2 : 2800 - 2000 A 

2000 - 1750 A 

1750 - 1350 A 

C0 2 : 2200 - 1750 A 

1750 - 1350 A 

0 3 : 7100 - 3000 A 

3000 - 2000 A 

2000 - 1350 A 


Photodissociation yield factors for oxygen 
and ozone — assumed to be unity 


Rate constants 

k 2 = 5.5 x 10 -34 (T/300) -2-9 cm 6 molec ~ 2 sec" 

-12 3 -1 

= 1.4 x 10 exp(-1500/T) cm molec sec 

k 5 = 3.0 x 10 -33 (T/300) -2-9 cm 6 molec " 2 sec ' 


Shardanand (ref. 18) 

Gast (ref. 19) 

Watanabe (ref . 20) 

Fabian & Libby (ref. 7) 
Watanabe et al . (ref. 21) 
Craig (ref. 4) 

Inn & Tanaka (ref. 22) 
Watanabe et al . (ref. 21) 


Kaufman 
(ref. 23) 
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ALTITUDE (km) 



NUMBER DENS 

Figure 1 - Comparative profiles of oxygen compounc 
the 5-reaction systems for a model atmc 
Mariner 4 occultation experiment. 



ALTITUDE (km) 
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